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ABSTRACT 

We present a grid of LTE atmospheric models and synthetic spectra that cover 
the spectral class range from mid-G to mid-K, and luminosity classes from V to 
III, that is dense in Tes sampling {AT^s = 62.5 K), for stars of solar metallicity 
and moderately metal poor scaled solar abundance ([^] = 0.0 and —0.5). All 
models have be e n comp uted with two choices of atomic line list: a) the "big" line 
lists of iKuruczl ( 1l992al ) that best rep roduce the broad-band s olar b lue and near 
UV fx level, and b) the "small" lists of lKurucz fc Peytremanru (119751 ) that provide 
the best fit to the high resolution solar blue and near-UV spectrum. We compare 
our model SEDs to a sample of stars carefull y selected from th e large catalog of 
uniformly re-calibrated spectrophotometry of lBurnashevI ( 119851 ) with the goal of 
determining how the quality of fit varies with stellar parameters, especially in the 
historically troublesome blue and near-UV bands. We confirm that our models 
computed with the "big" line list recover the derived T^s values of the PHOENIX 
NextGen grid, but find that the models computed with the "small" line list pro- 
vide greater internal self-consistency amo ng different spectral b a nds, a nd closer 
agreement with the empirical Teg scale of iRamirez fc MeleridezI (120051) . but not 
to the interferometrically derived Tog values of iBaines et all ( l2010l ). We find no 
evidence that the near UV ba nd discrepancy between mo dels and observations for 



Arcturus (a Boo) reported by lShort fc Hauschildtl ( l2003l ) and lShort fc Hauschildt 



( 120091 ) is pervasive, and that Arcturus may be peculiar in this regard. 



Subject headings: stars: atmospheres, fundamental parameters, late-type. 
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Introduction 



Critical comparisons of the observed absolute spectral energy distribution (SED), /a(A), 
to that predicted with PHOENIX computational models fror n the near UV to t he re d 
have been carried o u t for the S un, Procyon, and Ar c turus (see I Short fc Hauschildtl ( 120091 ) , 
Short fc Hauschildtl (120051 ) . and lShort &: Hauschildtl (120031 ) ). They presented evidence that 
both LTE and non-LTE models increasingly over-predict the fx level for A < 4000 A as 
the Teff value of the star decreases from that of early G stars to that of early K. The inves- 
tigations to date have been restricted to bright standard stars for which there is very high 
quality spectroscopic and spectrophotometric data, with the consequence that the quality of 
the SED fit has not been well sampled in the stellar parameter space T^q/ \ogg/[^]. Here, 
we take the first major step to rectify the situation by comparing a large grid of model SEDs 
spanning the cool side of the HR d iagram to observed SEDs taken from the spectrophoto- 
metric catalog of iBurnashevI ( 119851 ) . Our goal is to map out the goodness of fit, and the 
magnitude of any systematic discrepancies between model and observed SEDs as a function 
of the three stellar parameters, Tefj, log f?, and [^], with the more fundamental goal of con- 
straining the physical character of the sources of any discrepancies. We also compare our 
Tcff values inferred from SED fitting to other empirical and theoretical Tefj calibrations. 



2. Observed /a(A) distributions 



BurnashevI ( 119851 ) presented a large catalog (henceforth B85) of observed SEDs taken 



with photo-electric instruments on 0.5m class telescopes at various observatories in the former 
USSR from the late 19 60s to the mid 1980s and u niformly photometrically re-calibrated to 
the "Chilean system" (jShort fc Hauschildtl (120091 ) contains a more detailed description of 
the individual data sources included in this compilation). These data sets all generally cover 
the A range 3200 to 8000 A with AA = 25 A, and have a quoted "internal photometric 
accuracy" of ~ 3.5%. 



We have harvested from the B85 catalog a sample of stars that meet th e following criteria: 
1) ha ving an entry in the 5*^^ Revised Edition of the Bright Star Catalogue (JHoffleit &: Warren 
199ll ). (henceforth BSC5), 2) having a spectral class later than GO and earlier than MO, 3) 
having a luminosity class in the range from V to HI, 4) having an entry in the metallicity 
catalog of ICayrel et al\ (1200 ll ). 5) having no variability, chemical peculiarity, or binarity flags 
in BSC5, and 6) having an fx distribution that was not obviously inconsistent based on visual 
inspection, with those of other stars of the same BSC5 spectral class and within ±1 sub-class. 
The spectral and luminosity classes were taken from The Revised Catalog of MK Spectra 
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ypes for the Cooler Stars (JKeenan fc Newsoinll2000l ). the paper of iKeenen &: Barnbaum 
( 1999 ). or The Perkins Catalog of Revised MK Types for the Cooler Stars iKeenan &: McNeil 
( 119891 ) . in decreasing order of preference. For nine objects without spectral classes deter- 
mined by Keenan or his collaborators, we took spectr al types fro m the Genera l Catalogu e 



of Stellar Sp ectral Classificatio ns (Version 2010-Mar) ISkifl 



HD 147675 (Torres et al. 



222107 (IGrav et a/.ll2003h. H P 221673 (lAbtlll98lh . HD 156266, 115383 (IHarlan 



20061 ). (IHouk k Cowlev Il975l ). HD 71878 (IGrav et al. 



(J2010|): HD 50522 (Abtll2008l) 



2006) 



19741 ) 



HD 
HD 

34559, 49878 (I Adams et a/.lll935f ). For the six objects for which ISkiflfl (l2010f ) had multiple 
conflicting spectral classes, we took the class that agreed with that listed in BSC5. Appli- 
cation of criterion 6) necessitated only including stars for which there was more than one 
object of the same spectral class, to within ±0.5 sub-classes. We found that application of 
these six criteria, along with the metallicity criteria described below, limited us to 33 stars 
of spectral class ranging from G5 to K4 for giants and GO to G5 for dwarfs. All spectra were 
corrected for their heliocentric radial velocity, RV, using the RV values in BSC5. However, 
we expect the RV correction to have a very minor effect on the quality of spectral fltting at 
the low spectral resolution of the B85 data. 



We restricte d the sample to two metallicity ranges based on the value of [g] reported in 
the catalog of ICayrel et al\ (120011 ): a) stars of —0.1 < [^] < 0.1, and b) stars of —0.4 < 



111 



< —0.6. These two ranges are expected to contain stars of measured [^] approximately 



equal to 0.0 and -0.5, respectively, on the grounds that the quoted uncertainties, A[^], of 
determinations of oi^era// metallicity typically have quoted uncertainties of the order of ±0.1. 
For stars that contained multiple entries in the catalog of ICayrel et al\ ( 120011 ). we based the 



metallicity selection on either the mean, or the median, value of [jj], depending on whether 
there were obvious outlier values and skew. Metallicity range b) approximately represents a 
moderately metal poor population typical of the older thick disk, represented by stars such 
as Arcturus ([^]~ —0.7). Stars of [^]< —0.5 are increasingly rare in B85 catalog as [^] 
decreases, so -0.5 was the lowest [^] value for which we could harvest a significant number 
of stars of [jj] value within ±0.1 of each other. Applying the six criteria described above, 
we harvested 30 and 3 stars from the B85 catalog in each of the metallicity ranges a) and 
b), respectively. Table [1] contains the list of the st ars that were selecte d in each metallicity 
range, and the number of [^] measurements in the ICayrel et al\ (120011 ) catalog. Some stars 
had two independently measured SEDs in the B85 catalog, and where that is the case, both 
spectra are included in our analysis and provide a check on the internal consistency of the 
B85 spectrophotometry. The number of spectra is also indicated in Table [H 
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Figs. |2] and [3] show, for the hottest and coolest giant samples, the individual spectra of 
all stars of that sample, over-plotted with the sample mean spectrum and the spectra of 
the ±1(7 deviation from the mean. Note that the a values are most clearly meaningful for 
the two samples that contain more than a few spectra, namely, the G8 III and KO III stars 
of solar metallicity. For most samples, the spectra that passed our selection procedure fall 
within ±1(7 of the mean at most A values. We have separately calculated mean and deviation 
spectra for each sub-sample. 



2.1. KO III sample 

The KO III sample, shown in Fig. HJis special in that it contains more spectra than the 
other sub-classes (14 spectra of ten objects), and the distribution of spectra has a bifurcation 
that can be most clearly seen in the log A region from 3.52 to 3.55 (3300 to 3550 A). Therefore, 
we have broken up the KO III sample into two sub-samples, an / sub-sample of four stars (six 
spectra) with "low" , and an h sub-sample with six stars (eight spectra) with "high" , UV f\ 
level, respectively. The stars belonging to the / and h sub-samples are indicated in Table 
[U This bifurcation may be an artifact of our sorting of stars with an effective precision of 
about one sp ectral sub-class. Around K O, stars can be typed to half-sub-class precision (see, 
for example, iKeenan fc NewsomI (J2000[ )). and differences of ~ 0.5 spectral sub-classes would 
presumably have the greatest effect in the blue and near UV spectral bands. We note that 
two of the four stars in our "low" sample have a spectral class of K0.5 III. 



2.2. Arcturus 



To compare with the results of IShort fc Hauschildtl (120031 ) . it is necessary to include 
Arcturus (a Boo. HD 124897) (and other stars of the same SEP, wer e they available), in the 
current investigation. The metallicity catalog of ICayrel et al\ ( 1200 ll ) lists 17 measurements 
of [jj] for Arcturus ranging from -0.370 to -0.810, with most values in the range of -0.4 to -0.6, 
and the average being -0.54. Arcturus is the only star in the B85 catalog with a spectral class 
in the range from KO to K2 III with measured [^] values within ±0.2 dex of -0.5. Therefore, 
we cannot form a Kl III sample at [jj] Ri —0.5 for quality control inspection, or for forming 
a mean spectrum, as we have with the other spectral classes and metallicities. Moreover, 
Arcturus is flagged in BSC5 as being chemically peculiar, with the spectral class designation 
being K1.5IIIFe-0.5. 
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With an apparent V band magnitude of -0.04, the B85 spectrophotometric data for Arc- 
turus should be of relatively high quality. The B85 catalog contains three measured fx 
distributions, two of which are consistent with each other on the basis of visual inspection. 
The third is higher by about 0.05 dex for log A less than about 3.55 (A < 3550 A). We have 
formed a sample at spectral class K1.5 III for stars of nominal [^] equal to -0.5 using these 
three SEDs. 



3. Model grid 

3.1. Atmospheric structure calculations 

We have computed a grid of about 400 atmospheric models under the approximation 
of LTE with spherical ID geometry spanning a range in Tgfj from 6250 to 4000 K with a 
sampling, ATefj, of 125 K, and in \ogg from 1.5 to 5.0 with a sampling of 0.5, for [-^] values 
of 0.0 and -0.5. A AT^s interval of 125 K is close to the nomi nal Tpff difference between su c- 



cessive spectral sub-classes among GK stars (see, for example, iRamirez fc MelendezI (120051 )). 
Furthermore, we interpolate in our synthetic spectrum grid to achieve an effective AT^s of 
62.5 K, as described below. 



For spherical models, a value of the effective radius, Res, at ri2ooo = 1 is also required 
as input. However, we expect the value of Res will have a minor impact on the computed 
SED compared to T^s, \ogg, and [^], and we used values corresponding to a mass of IMsun 
at the given value of log g for all models. Because we are restricted to stars in the spectral 
class range G5 to K5, most of which are of luminosity class III (see Section [2]), we expect 
the mass distribution of the progenitor main sequence objects to be of the order of IMsun- 



Hauschildt et al\ (jl999bl ) found that for supergiant models of given T^s and log (7, which are 
even more spherically extended than our models, the choice of mass in the range 2.5 to 7.5 
Msun lead to only small differences in the relative SED. 

The observed SED will depend on the value adopted for the micro-turbulent velocity 
dispersion, ^x, which determines the extent to which spectral lines are blended in crowded 
spectral regions, such as the blue band, and therefore, affects the fx level to different extents 
in different A regions. We adopted values of 1.0 and 2 . km s~^ for logg values greater 



than and less than 3.5, respectively. IShort &: Hauschildtl (120051 ) found that ^t = 1 km s ^ 



provided a good fit to the solar SED from the near UV to the red. From analysis of line 
profiles in high resolution high S/N spectra of G and K III stars, previous investigators have 
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found ^T ranging from 1.5 to 2.0 km s~^ (lGravlll982l ). (JFovlllQTSl ). (JGustafsson et a/.lll974j ). 
In Fig. [T]we compare two synthetic SEDs computed with ^t = 1 and 2 km s~^ for a model of 
Tgff = 4875 K, logf7 = 2.0, and [^]=0.0, and a model of ^t = 2 km s~^ and the same stellar 
parameters except that T^s is 5000 K. The value of ^t makes an increasingly large difference 
as A decreases because of the increasing density of spectral absorption lines, and reaches a 
maximum of ~ 0.1 dex at 3200 A. This is a significant effect, and a proper exploration of 
the effect on the stellar parameters determined from SED fitting will require a model grid 
with a ^T dimension, and is beyond the scope of the current investigation. However, we can 
note qualitatively from Fig. [1] that models models in which ^t is over-estimated by 1 km 
s~^ will produce fits to the blue and near UV band flux level that over-estimate the value of 
Tcff by ^ 100 K. For all models, we adopted a mixing length parameter for the treatment of 
convective flux transport of one pressure scale height. 



For all models, we adopt scaled solar abundances. There is evidence that mildly metal 
poor stars of [^]= —0.5 have abundance distributions that are enhanced in a-elements with 



respect to the Sun by 0.2 to 0.3 dex (see, for example, iPeterson et al\ ( 1l993l )). However, 
given the scope of the model grid required for this initial investigation, we have decided to 
restrict ourselves to scaled solar ^ distributions. There has been recent c ontroversy over 



the va lues of the solar abundances for important elements such as CNO, with lAsplund et al. 



(I2OO4J ) and other papers in that Series, revising the values downward by 0.2 to 0.3 dex on 
the basis of 3p mod els that account for hydrodynamic turbulence. However, the work of 
Asplund et al\ (12004 ) has been found to be discrepant with abundances determined from 
very sensitive helioseismological fitting, and has recently been thrown into question on the 
grounds of self-consistency (see , for example, the very thorough recent investigation by 
Pinsonneault fc Delahayd (120091 )). Fo r simplicity, we restric t ourselves in this investigation 
to the solar abundance distribution of iGrevesse et al\ (119921 ). 



For each grid point, we compute two models that differ from each other in the choice 
of atomic line list. Our Series 1 and 2 models use the "big" and "small" line lists of 
Short fc Hauschildtl (120091 ) . respectively, and that paper contains a description of the con- 
tent of the two line lists, the evidence for and against each one, and an investigation of the 
effect that the choice of line list has on the computed SEP, esp ecially in the heavily line 
blanketed blue and near UV region. I Short fc Hauschildtl (120091 ) found that the choice of 
"big" or "small" line list has a significant impact on the computed fx values in the blue and 
near UV bands for both the Sun and Arcturus, and concluded that NLTE m odels of the Sun 
with t he "big" line list provide a better fit to the measured f\ distribution of iNeckel fc Labs 
(jl984j ). but that those with the "small" line list provide a better fit to three other measured 
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fx distributions, including one measured from space. Therefore, we have decided to evaluate 
the fit provided by models with both choices of line list here. We note that the choice of line 
list affects the model structure as well as directly affecting the synthetic SED calculation. 



3.2. Synthetic spectra 

For each model in both Series 1 and 2 we computed self-consistent sjTithetic spectra 
in the A range 3000 to 8000 A with a spectral resolution (-R = xx ~ 350 000) to ensure 
that spectral lines were adequately sampled. These were then degraded to match the low 
resolution measured fx distributions of B85 by convolution with a Gaussian kernel. B85 
present their data with a sampling of 25 A, but, neither B85, nor any of the original source 
publications that are still accessible, describe the instrumental spectral profile. Based on 
trial comparisons of observed SEDs with synthetic fx spectra convolved with Gaussians of 
various FWHM values, we found that a Gaussian instrumental profile with a FWHM value 
of 75 A provided the closest qualitative match to the spectral structure in the observed 
SEDs. Therefore, we have convolved all our synthetic spectra with a normalized Gaussian 
kernel of FWHM equal to 75 A. We note that this convolution also automatically accounts 
for macro-turbulence, which has been found to be around 5.0 km s~^ for G and K II stars 



Grayl (119821 ). We interpolate in log Teg between pairs of synthetic SEDs to obtain a grid with 
an effective sampling, ATefj , of 62.5 K. A Tes difference of 62.5 K at 3300 A corresponds to a 
difference in the log /a of a blackbody of 0.07, and this difference will be smaller for larger A 
and higher T^s values. For comparison, the discrepancy in n e ar UV band log fx level between 



LTE models and observations found by lShort fc Hauschildtl (120091 ) for Arcturus around 3300 
A is ?^ 0.15. 



All spectra, observed and synthetic, have been normalized by dividing by the average 
of their flux in the 6600 to 6900 A region, chosen to be just blue of the first signifi- 
cant telluric contamination bands due to O2 and H2O, to produce the distribution denoted 
/a,6750 = f (A=6750) " -^^8^. |5] and [6] show the comparison of the mean and ±la spectra of 
the observed /a,6750 distributions with the closest matching and bracketing synthetic /a,6750 
distributions for the hottest and coolest giant samples. Figs. [7] and M show the difference 
between the mean of the observed /a,675o distribution and the closest matching and bracket- 
ing synthetic distributions relative to the observed mean distribution, (/A,6750,MeanObserved — 
/A,6750,Modoi)//A,6750,MeanObserved for the Same two Samples. 



4. Goodness of fit statistics 

We interpolate all observed and synthetic spectra onto a uniform A grid of AA = 0.01 
A that moderately over-samples all spectra at all A values, and compute for each spectral 
class sample the root mean square relative deviation, a, of the mean observed /a,675o dis- 
tribution from the closest matching and bracketing synthetic /a,6750 distributions in the A 
range from 3200 to 7000 A, according to 



C" — "77 _^((/A,6750,Obs — /A,6750,Mod)//A,6750,Obs 



I 



where A^ is the number of A points in the A grid in the 3200 to 7000 A range. 

We also compute separate RMS values, o"biue and cXred, for the so-called "blue" and "red" 
sub-ranges of 3200 to 4600 A and 4600 to 7000 A, respectively. We also compute the mean 
relative deviations, Abiue and Aj-ed, for the two sub-ranges, according to. 



N 
I , 

A 



1 ^ 

-T7 /_^{fxfi750,Ohs — /A,6750,Mod)//A,6750,Obs (2) 

i 

A comparison of the blue and red a values indicates how well the synthetic spectra fit 
in the blue and near UV band given the quality of fit in the red band. The comparison of 
A between the bands allows an assessment of systematic discrepancies throughout the blue 
compared to the red because, unlike a, A retains its sign. A break-point of 4600 A was 
chosen on the basis of visual inspection of where the deviation of the synthetic from the 
observed spectrum starts to become rapidly larger as A decreases. In Tables [2] and [3] we 
present the a, cTbiuc, and cTred values for the Series 1 and 2 models, respectively, along with 
the best fit value of Tes and log g for each star. The value of the model [ jj] is also tabulated, 
although, its value was specified a priori rather than fitted. 



4.1. Trend with T^s 

Fig. [in] shows the variation of a, ablue, and (Jred with model T^s for the giant stars 
of solar metallicity. The best fit a value generally increases with increasing lateness of the 
spectral class. We note that the density of spectral lines generally increases with increasing 
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lateness. Therefore, this trend in the discrepancy between synthetic and observed SEDs could 
be explained by inadequacies is the input atomic data for bound-bound {b — b) transitions, 
or by inadequacies in the treatment of spectral line formation. The best fit a values for the 
Series 1 and 2 models differ negligibly for the earlier spectral classes, and there is marginal 
evidence that the Series 2 models provide a slightly better fit (lower a value) for the latest 
spectral classes. 



4.1.1. G5 III sample 

The behavior of the variation of (Jred with T^^ for the G5 III stars is peculiar and leads 
to a spurious result for the best fitted value of Tefj. From Fig. [7] it can be seen that this 
is caused by a broad absorption feature exhibited by the observed SED with respect to the 
model SEDs ranging from a log A value of 3.753 to 3.774 (5660 to 5940 A). As a result, 
the value of cXred is increased significantly, even for models that provide a good match to 
the overall spectrum. Therefore, our best fit value of Teg for the G5 III models is best 
determined from the blue band alone. This deficit of absorption in the synthetic SEDs with 
respect to the observed ones is consistently present in the individual observed spectra for the 
G5 III stars, spans 12 data points in the raw observed spectrum, and varies smoothly with 
wavelength over a range of 280 A. Therefore, it is likely caused by a cluster of spectral lines 
that are either missing, or are too weak, in the model spectra. We note that this discrepancy 
is either absent, or much less pronounced, in both the G4-5 V and G8 III stars, so appears 
to be localized in both T^s and log (7. 

We have examined histograms of the average numbers of spectral lines per A by atomic 
chemical species in a synthetic spectrum of a model of Tefj/log(y'/[^]=5250/2. 5/0.0, repre- 
sentative of the models that fit the observed G5 III SED. We compared three A ranges: the 
problematic 5660 to 5940 A range, and the bracketing ranges of 5000 to 5660 and 5940 
to 7000 A. All three ranges show approximately the same pattern and same absolute aver- 
age numbers of lines per unit wavelength for all atomic species that contribute a significant 
number of lines. We are unable to identify any suspect chemical species for which there is 
an excess or dearth of lines in the 5660 to 5940 region compared to neighboring regions that 
might provide a clue to the cause of the discrepancy. 

We have compared our mean observed spectra based on the B85 catalog for the G5 III and 
V, and G8 II I types to representative spectra from the stellar spectrophotometric library of 



Jacoby et al\ (119841 ). In Fig. [H] we show the comparison in the A range around the region 
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of the G5 III discrepancy. The iJacoby et all (Il984j ) hbrary does not have a spectrum for 
type G5 V, so we have compared our spectrum o f that class to spect ra for classes G4 and 
6 V. The G4 V star is designated 'TR A 14' and IJacoby et al\ (119841 ) flag it as a star that 
could not be identified. However, for lack of an alternative, we use it as a comp arison. There 



are sm all systematic differences between our mean B85 spectra and those of IJacoby et al. 



(Il984f) over a broad ran ge of A for all three spectral types. However, for the G5 HI sample. 



the iJacoby et al\ ( 1l984j ) spectrum is distinctly brighter than our mean B85 spectrum in the 
log A 3.753 to 3.774 range, as compared the the bracketing log A ranges, and is in greater 
agreement with our synthetic SED. We conclude that there may have been a data acquisition 
or calibration problem with the spectra in the B85 ca talog that is very pa rticular to the G5 
HI class. From Fig. [8] it can be seen that the B85 and iJacoby et al\ ( ll984j ) spectra are much 
more consistent with each other throughout this region for the G5 V and G8 HI classes. 
We note again that our B85 G5 HI sample suffers from small-number statistics in that it 
consists of three spectra of two stars. In what we follows we only draw conclusions for the 
blue spectral band of the G5 H I samp le. (A full systematic comparison of our B85 mean 
spectra with those IJacoby et al\ (Il984j ) throughout our A range is beyond the scope of this 
investigation, but, the comparison in Fig. [H]is generally encouraging.) 



4- 1-2. Red vs blue band 

For all samples, the best fit to the red band has a a value that is lower than that of the 
fit to the blue band by 0.05 to 0.1. This may partly reflect that all spectra were normalized 
to a common relative flux value in the red band (6750 A). However, from Figs. [7] and |9] it 
can be seen that the difference spectra show increasing variability around the zero line in 
addition to any systematic trend away from the zero line. This indicates that the quality 
of the fit to the mean observed spectrum worsens with decreasing A regardless of how the 
observed and synthetic spectra were normalized. For all solar metallicity giants, the Series 
1 models fitted to the blue band consistently give best fit Teg values that are one AT^s 
resolution element, 63K, higher than those fitted to the red. For the Series 2 models, both 
bands yield the same best fit Teg value for the KO to K2 HI stars, whereas for the two ends 
of the Teff range where both bands can be used, G8 and K3-4, the Series 2 models lead to 
the same pattern as the Series 1 modes, with the blue band fit yielding a higher Tes value. 
This may provide marginal evidence that the Series 2 models lead to a greater consistency 
of fit throughout the SED. 
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4.1.3. Series 1 vs Series 2 

For the samples at the two ends of our Tes range, G5 to G8 and K3-4, fits with the Series 
1 and 2 models both lead to the same value of inferred T^s, to within the ATgfj precision 
of the grid. However, for the intermediate samples, KO to K2, the Series 2 models fitted 
to the blue and to the overall SED consistently lead to fitted T^s values that are one ATgfj 
resolution element lower than that of the Series 1 models. Because the Series 2 models have 
less line blanketing, they predict greater flux in the blue with respect to the red than do 
the Series 1 models. Therefore, we expect them to yield lower fitted Tgg values to a given 
observed SED. 



4.2. Trend with \ogg 

Fig. [11] shows the fitting quantities for the G5 III, G4-5 V, and GO V stars of solar 
metallicity. Examination of these figures allows a limited assessment of how the quality of fit 
varies in the logg dimension at spectral class G5. The Series 1 models provide fits of similar 
quality to the total SED of G5 stars of both luminosity classes. III and V. By contrast, the 
Series 2 models provide a significantly worse fit to the class V than to the class III stars, with 
the quality of fit to class III differing negligibly from that provided by the Series 1 models. 
Furthermore, we note that this same discrepancy in the quality of fit provided by Series 1 
and 2 is also seen for the GO V stars. Fig. [11] shows that these trends in the quality of the 
fit to the total band is driven mainly by the quality of fit to the blue band. The suggestion 
is that the models computed with the "big" input line list provide a better fit to the blue 
band spectra of dwarfs than, and as good a fit to the giants as, the models computed with 
the "small" line list. 



For the G4-5 V stars, the Tes value fitted to the blue band with Series 1 models is 63 K 
higher than that fitted with Series 2 models, whereas, for the G 5 III stars, both Series yield 
the same blue band T^g value. The G4-5 V stars are anomalous among our solar metallicity 
stars in that the Series 2 models give a higher fitted T^s value to the red band than to the 
blue band. Unfortunately, our ability to compare the fit to the red band as a function of 
logg is undermined by the peculiarities with modeling the red band of the G5 III sample 
discussed earlier in this section. 
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4.3. Trend with [|] 

Fig. [T2] sliows tlie same quantities for tlie G8 III and Kl to K2 III stars of [^] equal to 
0.0 and -0.5, and allows a limited assessment of how the quality of fit varies in the metallicity 
dimension. For a given spectral class, both the Series 1 and 2 models provide a significantly 
better quality of fit to the stars of solar metallicity than to those that are metal poor. This 
may in part reflect the inappropriateness of scaled solar abundance models for fitting the 
population of stars of [g]= —0.5. By contrast, for a given T^s value, there is some evidence 
that all models provide a similar quality of fit, independent of [^]. We note that the metal- 
poor G8 III star has a best fit T^s value closer to that of the solar metallicity Kl to K2 stars 
than that of the solar metallicity G8 III star, and the value of a for the red and blue bands, 
and for the overaU SED, is similar for the 08 III/[|]= -0.5 and K2 III/[|]= 0.0 stars. This 
suggests that the strongest relation is the anti-correlation between T^s and a rather than 
that between [^] and a. 



4.4. Arcturus 

For the Arcturus sample (K1.5 III, [^]= —0.5) the Teg value fitted to the red band is 
higher than that fitted to the blue by 63 K. Therefore, models fitted to the red band would 
pr edict too much flux i n the blue and near UV bands . This is consistent with the results 



of IShort fc Hauschildtl (120031 ) and IShort fc Hauschildtl ( 120091 ) . who also compared models 



to the observed fx distribution of B85. We note that 62.5 K is the numerical precision of 
our TefT grid. Therefore, the most we can conclude is that the discrepancy in the best fit 
Tefj value between the red and blue bands is in the range of 31 to 125 K. By contrast, the 
solar metallicity giants of spectral type G8 to K4 do not show this trend, having best fit 
Teff values that are lower by 63 K in the red band (to within the precision of our model 
grid). Moreover, our only other metal poor sample, G8 III, also shows the opposite trend as 
Arcturus, with a red band T^q value that is 125 K (two ATefj increments) lower than that 
of the blue band. The suggestion is that Arcturus may be peculiar among giants generally 
in that models fitted to the red over-predict the blue band flux. 
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5. Comparison to other Teg calibrations 
5.1. Empirical calibration 



Ramirez fc MelendezI (120051 ) (RM05) present an empirical determination of T^g based 
on applying the Infrared Flux Method (IRFM) to 100 dwarfs and giants of spectral class 
from FO to K5 and [jj] from -4 to +0.4 as a function of many observables, inclu ding Johnson 



B — V color. This work is an extension of the definitive Teg calibration work of lAlonso et al. 



( )l999l ) and other papers in that series. They quote standard deviations of 30 to 120 K, and 
find that their IRFM T^q scale agrees to within 10 K with directly determined values for 
stars with measured angular diameters. RM05 do not present Tgfj as a function of spectral 
class. However, we have computed mean and RMS {a) B — V values for each of our spectral 
class samples using colors for individu al objects from the Catalogue of Homogeneous Means 



in the UBV System ( iMermillio dl 1 1 9 9 if ) . We then used the fitted Eqs., 1 and 2, and the fitting 
co-efficients of Tables 1 and 2, of RM05 to interpolate m B — V and produce empirical Teg 
values for each of our spectral class samples and metallicities. We checked our interpolated 
Teff values against, presumably, less accurate values found from linear interpolation in Tables 
4 and 5 of RM05 and found them to be consistent. In Table [Hand Figs. fT3] through [TSl we 
present a comparison of our T^q values fitted to our blue and red spectral ranges, and those 
of the RM05 calibration. Because of the peculiarities of fitting the red band of the G5 HI 
stars described in section HJ the red band Teg value has been suppressed in Figs. [13] through 



Baines et al\ (120101 ) (BIO) used the CHARA array to interferometrically measure the an- 
gular diameters of 25 K giants in the K band, including eight in the spectral class range KO 
to K4 HI with [jj] values within ±0.1 of 0.0, and two in the Kl to 2 HI range with [jj] values 
within ±0.1 of -0.5. They combined their angular diameter measurements with distances 
from Hipparcos and photometrically inferred bolometeric fiux values to derive T^s values. We 
have calculated averages for their Teff values for solar metallicity giants of spectral class Kl 
(three stars), K2 (two stars), K3 (1 star) and K4 (2 stars), and include these values in Table 
HJand Fig. [131 They also present a measured T^q value for an Arcturus analog (HD 170693, 
K1.5III-0.5). These results were pubUshed just as we were completing our investigation, so 
we included them for comparison. Unfortunately, none of the BIO stars are in the samples 
we selected from the B85 catalog. 
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5.1.1. Solar metallicity giants 

For the G8 to Kl stars, the models of both series fitted to the red band give Tes values 
closest to the RM05 values, and are consistently higher than RM05 by less than 100 K. For 
the K2 sample, our red band fit matches the RM05 value to within the ±62.5 K precision 
of our grid. The Series 1 models fitted to the blue band consistently yield Teg values that 
are higher than the RM05 Tcs calibration. An equivalent alternate interpretation is that the 
Series 1 models with the RM05 Tes value for their observed B — V color would predict too 
little blue band flux compared to the observed SED, while providing a closer match to the 
red band. 



The BIO Tcfj values for classes Kl and 2 III are is closest agreement with our Series 1 
values derived from the blue band, being an exact match at K2 III to within the precision of 
the respective values. Their value for the K3-4 stars is lower than our lowest value by ^ 50 
K. 



5.1.2. Solar metallicity dwarfs 

For the GO V stars, the Series 2 models fitted to the red band give a T^s value very close 
to the RM05 value. All models fitted to the blue band give Tcs values that are too high by 
p^ 100 K. For the G4-5 stars all models give Tcs values that are consistent with each other, 
but 100 to 150 K higher than that of RM05. These results are similar to those found for G5 
III stars in the blue band, so the discrepancy between our blue band Teff scale and that of 
RM05 at spectral class G5 seems to be independent of logg in the 2.5 to 4.5 range. 



5.1.3. Metal poor giants 

For the G8 III stars of [jj] = —0.5, the Series 1 and 2 models fitted to the blue band 
provide a very close match to the Teg value of RM05, better than the fit to the solar metallicity 
G8 giants. All models fitted to the red band provide T^s values that are 100 to 150 K lower 
than that of RM05. For Arcturus (K1.5III-0.5), all models fitted to the red band give T^g 
values in close agreement with RM05, being lower by 20 K. This result is consistent with 
our result for solar metallicity Kl and 2 stars, and suggests that the quality of our fit to the 
red band for early K giants is independent of metallicity in this [jj] range. The BIO value 
for their one K1.5 III star of [^] = —0.5 (HD 170693) is 75 K higher than our highest value. 
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5.2. Previous PHOENIX calibration 



Bertone et al\ ( 120041 ) (BBCR04) derived Teg values from fitting theoretical SEDs for 
models of [g] = 0.0 taken from the Nex tGen grid computed with an earlier version of 



PHOENIX (VIO) flHauschildt et a/.lll999bf ). The most important differences between VIO 
and V15 that are relevant to these models are 1) Improvements in the treatment of the line 
profiles so that many weak to moderate strength lines that were treated as Gaussian are now 
treate d with Voigt profiles, and 2) Improvement in the quality of a wide variety of EOS input 
data. (JHauschildt et a/.lll999al ) to SEDs of solar metallicity dwarfs and giant s of spectral class 
from early A to mid-M taken f rom the observational spectral libraries of iGunn fc Stryker 
( 119831 ) and IJacoby et al\ (Il984j ). We note that the NextGen models were computed using 
the "big" atomic line list of our Series 1 models. Among the main differences between the 
NextGen grid and the one presented here are: 1) the NextGen models of logg > 3.5 are 
computed with plane-parallel geometry, however, the effects of sphericity are expected to 
be negligible for that surface gravity range. 2) the NextGen models have the value of ^t 
set to 2 km s~^ throughout the grid, independent of logg. As discussed above, our value 
of ^T decreases to 1 km s~^ for the d w'arfs on the basis of our ex perience modeling the disk 
integrated flux spectrum of the Sun (IShort fc Hauschildtll2005l ). BBCR04 do not provide 
an interpolation formula for T^s, so we have interpolated linearly in their Table 1 to find 
Teff values for our spectral classes that are missing from their table. Table H] and Figs. [T3] 
through [IS] include a comparison the T^s values of BBCR04. 



5.2.1. Solar metallicity giants 



Generally, the BBCR04 Tes scale is 50 to 150 K hotter than the empirical RM05 scale 
in our spectral class range, and our models fitted to the red band consistently provide Tgff 
values that are lower than those of BBCR04 by ~ 50 K. We note that our Series 1 models 
fitted to the blue band provide consistently very good agreement to the NextGen Tes scale. 
That the Series 1 models would provide better agreement is expected since they are based 
on the same atomic line list as the NextGen grid, and the blue band is expected to show 
greater sensitivity to the treatment of line opacity than the red band. All of our fitted T^g 
values are lower than the BBCR04 value for the earliest spectral class, G5 III, which implies 
that the BBCR04 Tes scale would be even more discrepant with the empirical RM05 scale 
than ours is at the hot end of the sequence. 
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5.2.2. Dwarfs and metal poor stars 

Again, for the G4-5 V stars, our Series 1 models fitted to the hne-sensitive blue bands 
provide the closest match to the NextGen Teg value, which is what we expect. Somewhat 
surprisingly, for the GOV star all our models yield Tefj values that are 200 to 250 K lower 
than of the NextGen models. This result is surprising in that we expect discrepancies due to 
inconsistencies and inaccuracies in the treatment of line blanketing, among other modeling 
aspects, to be reduced with increasing Te^. We do not find such a gross discrepancy at GO 
V with the empirical RM05 T^e scale. 



6. Conclusions 

Our giant sequence samples T^^ well from spectral class G5 to K3-4. The a value of 
the best fit to the giant star spectra for all spectral bands increases with later spectral 
class, with the value for the overall SED increasing from ^ 0.05 at spectral class G8 III to 
^ 0.10 at spectral class K3-4 III. The quality of fit to the red band (4600 to 7000 A) is 
significantly better than that to the blue (3200 to 4600 A), with the best fit a value being ?^ 
0.02 for classes G8 to K2 III. (Unfortunately, a peculiar discrepancy between the observed 
and synthetic spectra unique to the G5 III sample prevents us from meaningfully including 
that sample in conclusions drawn for the red band.) Higher best fit a values are driven 
by an increasing variation around zero in spectra of the difference between observed and 
synthetic spectra as A decreases, and correlates with the increasing density of spectral lines 
in /a spectra with decreasing A. Therefore, the contrast between the best fit o"biue and cxred 
values probably reflects inadequacies in the input atomic line list and/or the treatment of 
line formation. 



For the Series 1 models, there is also a systematic effect with wavelength in that best fit 
models to the blue band yield a fitted Teg value that is one ATofj grid resolution element 
(63 K) higher than those fit to the red band. Equivalently, models fitted to the red band 
will predict too littl e flux in th e blue band. This re s ult is in stark contrast to the results of 



Short &: Hauschildtl (|2003[ ) and lShort &: Hauschildtl ( 120091 ) . who found that models fitted to 



the red band of Arcturus (K1.5 III) predicted too much blue band flux. We note that our 
Series 1 models fitted to the blue band provide Teg values that are within 50 K of the values 
derived with the PHOENIX NextGen grid, which also used the "big" atomic line list. 
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The Series 2 models ("small" line list) fitted to the blue band lead to best fit Tes values 
for the early K giants that are 63 K lower than do the Series 1 models ("large" line list), 
thus providing greater consistency of best fit Teg value among the spectral bands than do 
the Series 1 models. By contrast, for the two spectral classes for which we have reliable 
dwarf spectra, GO and G4-5, the Series 1 models provide a better fit to the blue band of 
the dwarfs than those of Series 2, while the two Series provide an equally good fit to the 
red band of dwarfs. We conclude that there is marginal evidence that the Series 2 models 
provide a better fit to solar metallicity giants while Series 1 models provide a better fit to 
solar metallicity dwarfs. 

When comparing giants of the same spectral class, the models provide a better fit at solar 
metallicity than at [^]= —0.5, However, when comparing giants of the same Tes value, the 
models give a similar quality of fit at both metallicities. Taking into account the metallicity 
dependence in the T^s calibration of the spectral classes, we conclude that the quality of fit 
provided by our models is not dependent, or only weakly dependent, on [^] in the range -0.5 
to 0.0. 



We find that Arcturus is peculiar among our stars in that it is the on/y object for which the 
best fit Teflf value to the blue band is lo wer than that for the red b and, by 63 K. This result 



i s qua litatively consistent with that of IShort &: Hauschildtl (120031 ) and IShort fc Hauschildt 



( 120091 ) ■ who found that models fit to the yellow and red bands of Arcturus yielded signifi- 
cantly more blue and near UV band flux than observed. Those authors concluded that there 
is an important source of near UV band continuous opacity missing from cool star atmo- 
spheric models. One goal of this work was to determine the extent to which this discrepancy 
was pervasive among late type stars. We conclude that it is not, and that Arcturus may be 
peculiar in this regard. If anything, there is a tendency for models fit to spectra of all other 
spectral classes to nnder-predict the blue fiux with respect to the red. 



For the solar metallicity giants, we find that the Series 2 models, which yield consistent Tgfj 
values between the red and blue bands, agree to within 100 K with the RM05 empirical T^q 
values for our full range of spectral classes. For the K2 to K4 III stars, we recover the RM05 
Toff values to within the AT^s resolution of our model grid. The agreement of the Series 2 
Tcff values with the RM05 values for the GO and G4-5 V stars is also very good. For the 
metal poor giants, the relative success of Series 1 and 2 models is mixed with G8 and K1.5 
III stars having contrasting results. By contrast, the very recently published interferometric 
Teff values of BIO are in very close agreement to the Series 1 models fit to the blue band for 
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early K giants. We conclude on balance that the Series 2 models ("small" atomic line list) 
provide greater internal self-consistency and agreement with empirical Tgfj scale. We advise 
that the large atomic line lists containing many theoretically predicted lines of Fe-group 
elements be used with caution. 



PH OENIX includes continuous molecular opacity from H2, Hg , H2 and MgH. lKurucz et al 



( 119871 ) found that including the photo-dissociation opacity of CH makes a detectable differ- 
ence to the predicted flux of the Sun in a localized region around 4000 A. Their calculations 
did not include the effect of line blanketing, so it is difficult to judge how detectable the dif- 
ference is. Nevertheless, because molecules become increasingly important as Tes decreases 
it wo uld be worthwhile investig ating the role of this opacity source in our models. In any 



case, 



Short fc Hauschildtl ( l2009l ) found that the continuous opacity in the Sun and Arcturus 
in the AA3000 to 4000 region is dominated, variously, by by H^ b — f, the combined b — f 
opacity due to metals, the combined f — f opacity of Mg and Si, and Thomson scattering, 
and the cross-sections for most of these are kown accurately enough that their dominance is 
not in doubt. 

CIS is grateful for NSERC Discovery Program grant 264515-07. The calculations were 
performed with the facilities of the Atlantic Computational Excellence Network (ACEnet). 
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Fig. 2. — G5 III sample (two stars, three spectra). Light gray dotted hnes: individual 
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Fig. 3. — Same as Fig. Ej but for the K3-4 III sample (four stars, four spectra) 
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and "red" bands. 
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Fig. 8. — Comparison of B85 and IJacoby et al\ (J1984J ) observed spectra. Solid black line: 
sample mean from the B85 catalog as in Fig. [21 solid g ray lines: syn t hetic spectra as in 
Figs. [5]and[6l dashed black lines: spectra from library of iJacoby et al\ (119841 ) for a G5 111 
star (BD+281885, upper panel), a G7 (HD 249240) and a G8 (HD 245389) 111 star (middle 
panel), and a G4 (TR A 14) and a G6 (HD 22193) V star (lower panel). Vertical lines: solid: 
region of discrepancy between the B85 and PHOENIX spectra for G5 111 stars; dashed: 
normalization region. 
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Fig. 9. — Same as Fig. [71 but for the K3-4 sample. Black lines: logg = 2.5, gray lines: 
log^f = 1.5. For clarity, the synthetic spectra the logg = 2.0 models have been omitted. 
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Fig. 10. — Giants of solar metallicity: Variation of a with model T^^. Solid line: Series 1 
models; dashed line: Series 2 models. Upper panel: Fit to total SED; Lower panel: Black 
lines: fit to blue band; gray lines: fit to red band. 
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Fig. 11. — Same as Fig. [TOl but for the dwarf stars. 
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Fig. 12. — Same as Fig. [TOl but for the metal poor giants. 
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Fig. 13. — Comparison of our best fit T^s values with other T^s cahbrations. Squares: Series 
1 models; Crosses: Series 2 models. Black symbols: Fit to the blue band; Gray symbols: fit 
to the red band. Dotted lines: Upper panel: Empirical calibration of RM05; Lower panel: 
PHOENIX NextGen models fitted to stellar spectral libraries (BBCR04). Triangles (lower 
panel): T^s values of BIO. 
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Fig. 14. — Same as Fig. [131 but for the dwarfs. 
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Fig. 15. — Same as Fig. [151 but for the metal poor stars. Note that BBCR04 only analyzed 
solar metallicity stars, so we only show a comparison to the Tes calibration of RM05. 
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Table 1. List of stars selected from B85 catalog for SED fitting. Spectral types are from 
the sources cited in Section [2], V values are f r om B SC5, and the [4] values are from 



Cavrel et all tOOlh 



Object 



Spectral 
type 



V 



Mean or median 



m 



Num Num 
[^] spectra 



HD50522 


G5 III-IV 


4.35 


0.05 


1 


1 


HD71369 


G5III 


3.36 


-0.05 


2 


2 


HD34559 


G8III 


4.94 


-0.10 


1 


1 


HD62345 


G8 Ilia 


3.57 


0.00 


1 


3 


HD100407 


G7III 


3.54 


-0.04 


1 


1 


HD115659 


G8 Ilia 


3.00 


-0.02 


2 


1 


HD 147675 


G8 KOIII 


3.89 


-0.05 


1 


1 


HD192947 


G8 Illb 


3.57 


-0.03 


2 


1 


HD71061 


K0.5 Illb 


4.51 


-0.04 


1 


2 


HD61935^ 


KOIII 


3.93 


-0.10 


2 


2 


HD62509^' 


KO Illb 


1.14 


-0.10 


7 


2 


HD123139^ 


KO Illb 


2.06 


0.04 


2 




HD1253511 


KOIII 


4.81 


-0.08 


2 




HD177241^ 


G9 Illb 


3.77 


0.03 


2 




HD188947^ 


KOIII 


3.89 


0.05 


2 




HD2055121 


K0.5 III 


4.90 


0.05 


4 




HD216228^ 


KOIII 


3.52 


0.03 


3 




HD219449I 


KOIII 


4.21 


-0.08 


2 




HD62044 


KlIII 


4.28 


-0.02 


1 




HD71878 


Kl III 


3.77 


-0.01 


1 




HD206952 


Kl III 


4.56 


0.04 


1 




HD156266 


K2III 


4.73 


-0.03 


1 




HD161096 


K2III 


2.77 


0.06 


2 




HD49878 


K4III 


4.55 


0.05 


1 




HD81817 


K3 Ilia 


4.29 


0.09 


1 




HD211073 


K3III 


4.49 


-0.07 


2 




HD221673 


K4 Illb 


4.98 


-0.03 


1 




HD115383 


GOV 


5.22 


0.08 


5 




HD141004 


GOV 


4.43 


-0.02 


6 




HD20630 


G5 V 


4.83 


0.01 


3 
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Table 1 — Continued 



Object 


Spectral 
type 


Mean or median 

V \H 


Num 

rAi 


Num 

spectra 


HD117176 


G4 V 


4.98 -.08 


4 


1 



HD138905 G8.5 III 3.91 

HD222107 G8 III-IV 3.82 

HD124897'^ K1.5 IIIFe-0.5 -0.04 



0.41 


2 


1 


0.49 


2 


2 


0.54 


17 


3 



iKO III stars of "low" (1) UV f\ level; see text. 
^KO III stars of "high" (h) UV fx level; see text. 



^Arcturus, a Boo 
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Table 2. Series 1 models: Closest match models to mean sample spectra and goodness of 

fit statistics. 



Spectral type 



Total SED Blue 

Teff (a) logg Teff (a) 



Red 



m 



G5III 




a 


2.5 


5312 


;0.058) 




a 


0.0 


G8III 


5125 


;0.039) 


2.5 


5125 


;0.058) 


5062 


;o.oi3) 


0.0 


KOIII 


4875 


;0.049) 


2.0 


4875 


;0.077) 


4812 


;o.oi2) 


0.0 


KO III^ 


4938 


;0.046) 


2.0 


4938 


;0.069) 


4875 


;o.oi6) 


0.0 


KO IIJi 


4812 


;0.059) 


2.0 


4812 


;0.093) 


4750 


;o.oi6) 


0.0 


Kl III 


4688 


;0.079) 


2.5 


4688 


;0.125) 


4625 


;o.oi4) 


0.0 


K2III 


4562 


;0.062) 


2.0 


4562 


;0.098) 


4500 


;o.oi7) 


0.0 


K3-4 III 


4250 


;o.io9) 


2.0 


4250 


;0.169) 


4188 


;o.o4i) 


0.0 


GOV 


5938 


;0.039) 


4.5 


5938 


;0.059) 


5812 


;o.oi4) 


0.0 


G4-5 V 


5688 


;0.049) 


4.5 


5688 


;0.075) 


5625 


;o.o2i) 


0.0 


G8III 


4688 


;0.063) 


2.5 


4688 


;0.094) 


4562 


;o.oi3) 


-0.5 


K1.5 III^ 


4250 


;0.095) 


2.0 


4250 


;0.147) 


4312. 


;0.026) 


-0.5 



^See text. 

^Arcturus, a Boo 

^KO III sub-sample of "high" UV fiux (see text). 

^KO III sub-sample of "low" UV fiux (see text). 
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Table 3. Series 2 models: Same as Tabled 



Spectral type 



Total SED Blue 

Teff (a) \ogg Teff (a) 



Red 

Teff (cr) 



m 



G5III 




a 


2.5 


5312 


;0.065) 




a 


0.0 


G8III 


5125 


;0.037) 


2.0 


5125 


^0.057) 


5062 


;o.oi2) 


0.0 


KOIII 


4812 


;o.o5i) 


2.0 


4812 


;0.082) 


4812 


;o.oi2) 


0.0 


KO III^ 


4875 


;0.046) 


2.0 


4875 


;0.074) 


4875 


;o.oi4) 


0.0 


KO III' 


4750 


;0.057) 


2.0 


4750 


;o.o9o) 


4750 


;o.oi8) 


0.0 


KlIII 


4625 


;0.074) 


2.5 


4625 


;o.i2i) 


4625 


;o.oi5) 


0.0 


K2III 


4500 


;0.059) 


2.0 


4500 


;0.095) 


4500 


;o.oi8) 


0.0 


K3-4 III 


4250 


;o.ioi) 


1.5 


4250 


;0.155) 


4188 


;0.036) 


0.0 


GOV 


5938 


;0.048) 


4.5 


5938 


;0.076) 


5875 


;o.oi2) 


0.0 


G4-5 V 


5625 


;o.o6o) 


4.5 


5625 


;0.094) 


5688 


;0.022) 


0.0 


G8III 


4688 


;0.058) 


2.5 


4688 


;0.086) 


4562 


;o.oi3) 


-0.5 


K1.5 III^ 


4250 


;0.096) 


2.0 


4250 


;0.149) 


4312 


;0.027) 


-0.5 



^See text. 

^Arcturus, a Boo 

^KO III sub-sample of "high" UV flux (see text). 

^KO III sub-sample of "low" UV flux (see text). 
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Table 4. Comparison with empirical T^s calibrations of RM05 and BBCR04. 



Spectral type B-V (a) 



Series 1 Series 2 

Blue Red Blue Red 



RM05 BBCR04 BIO 



G5III 


0.882 


'0.019 


) 5312 




5312 




5137 


5410 




G8III 


0.930 


'0.004 


) 5125 


5062 


5125 


5062 


4964 


5123^ 




KOIII 


1.043 


'0.002 


) 4875 


4812 


4812 


4812 


4721 


4870 




KO III^ 


1.018 


'0.011 


) 4938 


4875 


4875 


4875 


4781 


4870 




KO III^ 


1.080 


'0.000^ 


4812 


4750 


4750 


4750 


4650 


4870 




Kl III 


1.115 


'0.003^ 


4688 


4625 


4625 


4625 


4592 


4710^ 


4737 


K2 III 


1.160 


'0.006^ 


4562 


4500 


4500 


4500 


4531 


4550 


4562 


K3-4 III 


1.408 


;o.oi4; 


4250 


4188 


4250 


4188 


4118 


4243^ 


4134 


GOV 


0.595 


'0.004 


) 5983 


5812 


5938 


5875 


5864 


6140 




G4-5 V 


0.695 


;o.oio 


) 5688 


5625 


5625 


5688 


5519 


5670 




G8 III-0.5 


1.010 


'0.000 


) 4688 


4562 


4688 


4562 


4684 






K1.5 111-0.5^ 


1.211 


;o.oo9; 


4250 


4312 


4250 


4312 


4332 




4386 



^Tcs values found from linear interpolation in Table 1 of BBCR04. 

^Arcturus, a Boo 

^KO III sub-sample of "high" UV flux (see text). 

^KO III sub-sample of "low" UV flux (see text). 



